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INTRODUCTION 


This report summarizes accomplishments of a Factory Mutual 
Research Corporation (FMRC) project on the Prediction of Fire 
Dynamics for the NIST grant period August 1991 through August 1992. 
Work performed under a subcontract by Professor H.W. Emmons on 
Ceiling-Jet Dynamics and on Development of Strategies for Perform-— 
ance Fire Codes is first described under Task 1. The accomplish- 
ment of three tasks performed at FMRC are then presented in 
summaries of Tasks 2 — 4. All of this work is aimed at the devel- 
opment of subroutines or algorithms that can be used in NIST/BFRL 
comprehensive computer models. | 

During the past year, there has been further progress in the 
development of predictive models for flame radiant heat flux and in 
the development of a practical laboratory test method for measuring 
the smoke point of solid noncharring or charring materials. This 
progress should allow continuing improvements in the accuracy and 
applicability of fire propagation theories. 


Task 1. Prediction of Fire in Buildings (H.W. Emmons, Gordon McKay Professor 
of Mechanical Engineering Emeritus, Abbott and James Lawrence Professor 
Engineering Emeritus) 

The objective of this task is to advance the knowledge of the science of 
fire physics, chemistry, psychology, and physiology, so as to make rational 
quantitative design for fire safety possible. 

The two major items worked upon this year were: 

The simplified theory of ceiling jets. 


2. The advancement and use of fire models in fire safety engineering. 


1. The Ceiling Jet 

Present knowledge of the flow of hot fire gases along a ceiling consists 
of two parts. Practical engineering generally uses empirical data from some 
source judged to be equivalent or a field model, used in a few cases where 
nothing else was possible. The three dimensional field model would be expect- 
ed to agree excellently with experiment but I know of no detailed check. It 
is essential for research but too slow and too costly for most practical work. 

The wide spread use of field models depends upon future computer cost and 
availability. They will certainly not be used if there are sufficiently ac- 
curate and simpler methods available. We should note that the flow through a 
nozzle is never computed by a field model. A simple application of Bernoul- 
li's equation plus an empirical coefficient correcting for boundary layers, 
turbulence, and detailed nozzle design serves all practical purposes. This 
Task #1 was undertaken to provide a simple method for the ceiling jet. 

The conservation of mass, momentum, and energy with assumed top hat velo- 
city and temperature profiles provide what would be expected to be an adequate 
theory. The success of this expectation ws presented at the third IAFSS meet- 
ing at Edinburgh (1) (and Home Fire Project Report #82). For friction and 
entrainment the theory is adequate (although it cannot always meet a general 
inlet nor outlet condition) just as is true for hydraulic channel flow. For 
heat transfer, the solution cannot meet required boundary conditions on either 
end and does not agree with available experimental results. The inlet condi- 
tions are perhaps not too important, but during the transient phase and for an 
open ended corridor or ceiling the Richardson number must be 1. This has been 


a most difficult problem which only recently has shown some hope of solution. 


If we look at the simplified x momentum equation, it is not obvious that 
density increase by heat transfer should cause so much trouble. The published 
reports (1) contain a somewhat complex solution. A far simpler one has re- 


cently been found as follows: 


2 
mu sa =E(py - p) 66, + aS Ry (Ref. 1 equation 2) 1 


The right side can be simplified to 
. A, -£ 2 
mu, = 6 (8° (,-0)), 2 


which now integrates immediately to 


{~ 2. =E (o,-0) 6° + Cc 3 
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Evaluating the constant of integration at the position where the flow is crit- 


ical (subscript c) i.e., where Ri = 1 


. 3u 
e235 Mm 
and equation 3 reduces to the simple relation 
zy pl 
Rie ta ate 2 
Cc 


This relation appears hopeful for the existence of a tranquil flow upstream of 
Ri = 1, because a decrease of velocity from the critical value u,/u > 1 
implies an increase of Ri about 1 which is desired. However if u and u, are 


eliminated from equation 5 in terms of Ri, we get 


3 T - dba 


Te feted Mie 
Since the left side has a maximum at Ri = 1 (see Figure 3 of #82 of Ref 1), a 
real solution of Ri exists only for T < T, i.e. locations downstream of Ri = 
ae 
To see what is happening, equation 3 is expressed in terms of n = 6/6, 


Ra = ®a/Par R = 0/0, 
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This has been plotted in Figure 1 for three cases 
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We see that upstream of Ri = 1 there is only one real root and this at nega- 
tive depth. Downstream of Ri = 1 there are three real roots one with tranquil 
flow, one with shooting flow and an unreal negative depth. Within the top 
hat theory there are no solutions which can describe a ceiling jet with 
assorted downstream conditions. 

To see how nature, in the form of Chobotov's careful experiments (2), 
manages to avoid the large increase of Ri predicted by the simple theory, the 
terms of the x momentum and energy equations (8,9) were evaluated by plotting 
yee eh ean) agand! (Tt= 1.) 


: pa 2 2 
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Me, (en(T - ed) = -h 9 


and graphically measuring their slopes and thus evaluating the effective f and 
graphically measuring their slopes and thus evaluating the effective f = 
1/1/2pu° and SU: = q /ouc, aT. The results are given (£n-&n) as functions of 
Reynolds number on Figure 2. The values are nearly constant (f = .0652, St = 
-0408) for the first two meters of the experimental 5 meter channel and then 
fell to f = .046, St = .030. The values of the friction factor are much larg- 
er than seems reasonable for these boundary layer flows. These high values 
probably include some of the effects of the top hat assumption and the fact 
that the experiment did satisfy all the boundary conditions which the simple 
theory cannot do. 

Where things stand now is that practical work must either continue to com- 
pute the ceiling jet by a field theory or use the available empirical results 


from some source judged to be equivalent. 


AN ENHANCED TOP HAT THEORY 

One might think that the next more complex theory to the simple top hat 
assumption might be one with a more realistic velocity and temperature pro- 
files. These changes cannot improve the ability to meet inlet and outlet or 
leading edge requirements. This was inadvertently shown by Chobotov in his 
thesis not to be true. He developed such a theory which failed to agree with 
his experiments. To satisfy the required boundary conditions it is essential 
for the momentum equation normal to the ceiling y to force required depth 
changes to meet inlet and outlet conditions. 

Although a few of the resultant formulas have been given in progress 
reports, no hint was given as to how they were obtained. This will be briefly 
done in the following. 

For the present we assume steady flow and neglect the effect of friction 


on the y flow component. The basic momentum equations are: 


(ou) , - (euv)y =Pyeare 10 


2 
(pv dy + CSV) “Py + 06 11 
The top hat assumptions made are: 
1. u, T independent of y to depth 6 
2. v= uéd, y/é 12 


which assumes v to be a component of u produced by the change of depth of the 
layer surface which then decreases linearly to zero at the ceiling. Substi- 


tuting 12 into 11 gives for the vertical pressure gradient within the ceiling 


jet, 
2 2 2 
OU Grek) ree pu 6 y 
XLV X 
py = og - [ x einer 13 


of which the first term on the right is the usual hydrostatic effect while the 
next two terms are y momentum adjustments. The pressure throughout the ceil- 
ing jet is found by integration of equation 13 from y to 6 (since the pressure 


is known at 6 but not at the ceiling.) 


2 
-p(x,y) = -p(é6) + pg(6 - y) - (pu~s8,), SS ort LTO 14 


note that the last two terms of equation 13 have been combined after the y 
integration. The x momentum equation now becomes 
aii é y 2 2 
2 “ z ic Comme Va) 5 
Cour ye (= (ey) Fagleey) yo) (Couns oe ), °F 15 
and now the x momentum equation for the ceiling jet is obtained by integrating 
equation 15 in y from 0 to 6(x). After some effort to combine terms the final 


xX momentum equation for the ceiling jet is 


2 ave je 
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2 
(pus), = 5 3 XX 5 16 
This form of the momentum equation is simpler than previously published and 
has not yet been extensively explored. 

For the case of no friction, heat transfer, nor entrainment, equation 16 


can be integrated. The differential equation becomes: 


2 
RASL NG) uki2 0 17 


1 
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where n = 6/6,, Z = x/&,, and subscript c refers to the critical condition. 


Equation 17 has the solution 
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and the integration starts at z = L/6é.; the open end of a corridor or the 


a? 
leading edge of a transient ceiling jet. 

If the integration is started at the end and computes the ceiling jet 
depth back to the source which is far away, the denominator must approach 
zero. Since the denominator contains a square root, a zero denominator is not 


sufficient since the integration of a reciprocal square root is finite for all 


finite ye sthuseco have a ceiling jet of small finite depth back to a distant 
source requires that the denominator and its derivative must both be zero. 


Thus for a finite source of depth n1 


Cy = 3/nl + 1.5 n1° 21 


Cp = -cynl - 3an nl + .5n19 22 


The values of cy, Cp and the ceiling jet slope and curvature at the critical 
point are plotted in Figure 3. The most important result is the value limits 


set for the critical ceiling jet slope (ng) There is a real slope for 


ee 
shooting flows, ni < an but only imaginary values for tranquil flow. This is 
the effect of the square root. This limitation of the flow to shooting flows 
is in agreement with a remark by Zukoski who said that the Richardson number 
behind the front of an advancing ceiling jet was usually slightly less than 
one. 

If we now look at the source of a ceiling jet flow we are given in addi- 
tion to the mass flow rate and density, the depth and the slope of the surface 
(+5) L- The x momentum equation can be solved exactly as before resulting in 
the solution: 


x Lae ES On Rwers yt Gtete aa: 
oe J a Ve 23 


ie Fe n 
inate (c, + Cyn -32n yn - 3) 


Where c, and cy are again set by the long range ceiling jet depth, nl, as in 


equations 21 and 22. However they are not related to the inlet conditions. 


C4 = Ciccone aa 3/n; + 1.5 (Hg ay 


C5 ipva?e Teo ae Gi = cny 25 
Again to have the solution contain a finite depth at long range requires that 
the values of c, and cy be the same as before. Thus we can see the complete 
solution process as follows. 

The boundary conditions to be satisfied are; critical flow at z = L/é,; 
and?7 t= Nis Qe = (no)y at z2 = 0. By inserting equation 21 and 22 into 25 we 
get 


ee 


2 n1 me) 


il To. = 50UN ni/n; epee ey gs es + 3n1°) (n, = %y1) 26 


from which n1 is computed. Cc, and cy are then found from 21 and 22. These 


constants are then used to compute upstream from z=L/é, by equation 18 and 


downstream from z = O by equation 23. The two solutions are then connected 


together somewhere in between. 


2. Past Activities and Publications Related to the Advancement and Use of 


Fire Models in Fire Safety Engineering (See Item 7 for an explanation of 


the use of Table 1 in fire models) 


The conference on Fire Safety Design in the 2ist Century May 8-10, 
1991 at Worcester Polytechnic Institute. Presented the paper 
"Strategies for Performance Codes in the U.S." 

Third Symposium of the IAFSS at Edinburgh July 8-12, 1991 Presented 
the paper "The Ceiling Jet in Fires". 

Participated in the NFPA Model Users Group meeting October 2, 1991 at 
Worcester Polytechnic Institute. 

Presented a seminar on "Fire Science and the New Fire Safety Engin- 
eering" October 16, 1991 at Chemical Engineering Seminar Series at 
MIT. 

Participated in the Annual Fire Research Conference October 28-30, 
1991 at NIST. 

Participated in the CIB meeting October 30 - November 1, 1991 at 
NIST. Presented the paper “Validating Models and Algorithms". 
Presented a tutorial on “Fire Model Applications - the MGM Fire" 
January 15, 1992 at NFPA Model Users Group meeting at NFPA. Since 
this was a tutorial, it was important to discuss how to find the 
model input data and the manner of defining burning objects as clus- 
ters of numerous items. A copy of the available input data constants 
distributed at the meeting is attached as Table 1. 

Presented a seminar "The New Fire Safety Engineering" March 26, 1992 
at the Mechanical Engineering Seminar Series and the Center for 


Applied Energy Research at the University of Kentucky. 


9. Presented a lecture on "Genesis and Future of Fire Modeling" April 4, 
1992 as the first lecturer on the new Howard W. Emmons Distinguished 
Lecture at Worcester Polytechnic Institute. 

10. Presented a seminar “Fire Science and the New Fire Protection Engin- 
eering" April 29, 1992 at the Thermal Sciences Seminar Series 


Mechanical Engineering Department MIT. 


REFERENCES 

1. H. Emmons, "The Ceiling Jet in Fires," Third International Symposium on 
Fire Safety Science - also the more detailed version Home Fire Project 
Technology Report #82, 1991. 


2. Chobotov, M.V., Gravity Currents with Heat Transfer Effects, Ph.D. Thesis 
California Institute. of Technology, 1987. 
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acetone 
CH3-C-CH3 
" 

) 
re) density kg.n3 813 
k thermal conductivity w/m°K .19 
Cy specific heat j/kg°K 2114 
Ho #107! combustion heat j/kg 3.08 
H*10-© pyrolysis heat j/kg - 
Tig Ignition temperature°K 834 
Tp Flame temperature°K 1400 
ke air fuel ratio (stoi) 9.45 
Y air fuel ratio 16.6 
€ emissivity 4 
K absorption coef. m7! 1 
X fraction heat release 91 
XR Fraction heat radiated ES 
A burn rate —— 
f(CO5) fraction pyro C05 1.734 

n 
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f(HC) FRACTION PYRO (HC) de 
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F(H50) fraction pyro H50 al 
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Task 2. Models for Turbulent Flame Chemistry and Radiation 
COLO MR Use il) OL OLL) 


Radiation from turbulent diffusion flames controls the burning 
rates of moderate and large-scale fires (1) and is, therefore, of 
central importance to the flammability problem. Recently, we have 
developed a global flame radiation model (2) which predicts the 
radiant fraction from buoyant turbulent (fuel-jet) diffusion flames 
in terms of the fuel’s laminar flame smoke-point value. The model 
considers the separate roles of soot and gaseous radiation. The 
effective flame radiation temperature is predicted by the model in 
terms of the: 1) incompleteness of combustion, 2) radiant heat loss 
and 3) turbulent mixing, resulting in good agreement with experi- 
ment (3). 

It is argued that the soot absorption coefficient (or soot 
volume fraction) is: 1) proportional to the turbulent macro-scale 
flow time in the soot formation region of the flames, and 2) 
inversely proportional to the chemical time for soot formation. 

The chemical time for soot formation is shown to be approximately 
proportional to the smoke-point heat release rate on the basis of 
the laminar flame smoke-point models of Kent (4) and Guilder (5). 
Markstein’s general correlations (6,7) of smoke-points in terms of 
the stoichiometric oxidant to fuel mass ratiof, S, and adiabatic 
stoichiometric flame temperature, T,,; for nitrogen diluted fuel and 
oxidant mixtures show that: 1) the soot formation is quite sensi- 
tive to S which controls the local velocity of the gas flowing 
through the diffusion flame toward the fuel side, while 2) the soot 
formation is relatively insensitive to T,,. The latter result 
allows us to generalize the turbulent flame radiation model to 
include nitrogen diluted fuel mixtures and nitrogen diluted/oxygen 
enriched oxidants. 

The predicted radiant fractions are in excellent agreement 
with experiment for different oxidant and fuel mixtures covering a 
wide range of flame sootiness. The model is tested against both 
axisymmetric (8) and planar (9) free-jet buoyant turbulent diffu- 


sion flames over a range of heat release rates. We are currently 


16 


applying the model to wall fire flames which have smaller radiant 
fractions due to: 1) convective heat loss to the surface, 2) 
significant dilution of the supplied fuel gases by products of 
combustion diffusing back toward the fuel source, and possibly 3) 
reduced macro-scale mixing times due to wall shear stress. 

The model requires knowledge of the fuel smoke point which 
effectively characterizes the fuel’s chemical structure. Methods 
are now available for evaluating the required smoke point of 
gaseous, liquid, and solid (charring and non-charring) fuels of 
interest. 

By developing algebraic curve-fit formulas for the soot and 
gaseous emissivities, the model is self-contained and can be 


readily implemented using spread-sheet software. 
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Task 3. Turbulent Flame Heat Transfer to Surfaces (G.H. Markstein) 
The scientific understanding of wall-fire combustion and in 


particular the development of pertinent theoretical models (1-4) 
require experimental data on flame radiation so that these models 
can be applied to fires of hazardous scale where radiative energy 
transfer plays a controlling role. The purpose of the present work 
is to extend the range of such data beyond the upper limits of fuel 
sooting tendency and flame height attained in previous investiga-— 
tions (2,5-8). In the approach used here, wall fires of solid 
fuels are simulated by burning gaseous hydrocarbon fuels under 
steady-state conditions. The four fuels, methane, ethane, ethylene 
and propylene, are used to cover a sufficiently wide range of 
sooting tendency. 

Initially, the study concentrated on the simulated overfire 
region by releasing the entire fuel flow from a horizontal slot 
burner placed adjacent to a vertical water-cooled metal plate 
(9,10) while more recent work (11) is concerned primarily with the 
pyrolysis region, simulated by supplying gaseous fuels from a 
water-cooled vertical sintered-porous-metal surface. The burner 
developed for this study consists of ten 132-mm-high and 380-mm-— 
wide vertical porous-metal panels topped by a 660 mm high water- 
cooled solid-metal heat transfer plate divided into five similar 
segments. To approach two-dimensional flow conditions as closely 
as possible, water-cooled side walls of 150 mm depth are attached 
to the burner over the entire height of 1980 mm. 

The water-cooling passages embedded in the porous metal and 
the heat transfer plate are instrumented with differential 
thermocouples for measurement of the total heat transferred to each 
panel and segment by the flame. Other instrumentation includes a 
wide-view-angle radiometer measuring total radiant emission and a 
slit scanning radiometer to determine the vertical distribution of 
radiant emission by horizontal slices across the wall fire. 

In the slot-—burner work (9,10) it was found that the radiative 
fraction of total heat release rate was reduced significantly by 
placing the burner adjacent to a water-cooled wall, compared to 


18 


values obtained with free-burning slot-burner or jet flames. This 
work also yielded dimensionless correlations of radiant emission as 
a function of height independent of heat release rate, by introduc- 
ing a radiation flame-length parameter. The latter varied with the 
1/2 power of heat-release rate, suggesting that flame heights are 
controlled primarily by radiant cooling rather than turbulent 
mixing, and implying a mean volumetric heat-release rate indepen- 
dent of fire size and fuel type. 

Work with the porous-metal burner (11) showed that for any 
fixed fuel mass flux, the radiance N, increases linearly with 
height, z, after an initial jump at the flame base. At fixed z, 
linear relationships between reciprocal radiance N,’ and reciprocal 
mass flux m"'} were obtained. For methane and ethane these rela- 
tionships held over the entire range of mass flux, while, for 
ethylene and propylene, an abrupt transition separating two regimes 
occurred at critical mass fluxes. The ordinate intercept N., 
representing an asymptotic radiance, and the slope s of the linear 
relationship, are functions of height z. WN. and s?* increase with 
fuel sooting tendency at fixed z. 

Heat fluxes to the porous-metal panels decreased with increas-— 
ing mass flux near the flame base, but were nearly independent of 
mass flux near the top portion of the porous-metal burner. These 
heat fluxes increased with fuel sooting tendency except for a 
reversed trend between methane and ethane. At sufficiently high 
mass flux, a jump of heat flux occurred at the junction between the 
porous-metal burner and the solid-metal heat transfer plate, and 
the heat flux increased further with height in this over-fire 
region. In this zone, the data were separated into a lower over- 
lapping set for alkane fuels, methane and ethane, and a higher 
overlapping set for the alkenes, ethylene and propylene. 

For a more complete understanding of wall-fire heat transfer, 
additional data are needed. In particular, separation of total 
heat flux to the wall into its radiative and convective components 
is an essential requirement for the development of a predictive 


theoretical model. Instrumentation for performing pertinent 
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measurements is currently being developed. However, preliminary 
estimates of the separation can be obtained from the current data. 
Figure 1 shows a plot of qt vs radiance N, at a fixed mass flux m" 
and a fixed height z for the four fuels. Similar linear plots were 
obtained for other values of m" and z. The ordinate intercept, 
which represents the heat flux extrapolated to zero flame radiance, 
may be interpreted as the convective component of heat flux, while 
the slope multiplied by the radiance may be regarded as the radia-— 
tive component for each fuel. 

Intercepts as a function of mass flux for two heights are 
shown in Fig. 2 and Figure 3, which show that they decrease with 
increasing mass flux. This decrease is more pronounced at a height 
of 198 mm than at the greater height of 726 mm. A theoretical 
argument’* shows that the slope should vary from a value of 1.8 f& 
for optically thin to a value of m for optically thick isothermal 
flames. Figure 4 shows that at a height of 726 mm, the slope 
varies from 1.6 & (near optically thin) for low mass flux of 5 
g/m*s ‘to “about 1.0 1% (optically thick) for highsmass fluxivotes 
g/m*s, while Fig. 5 shows that at a mass flux of 15 g/m’s, the slope 
maintains the optically-thick value of 1.0 ® over the entire height 
range. The new instrumentation will permit verification of this 
tentative interpretation, and may, in addition, indicate possible 
deviations from isothermal conditions, as well as blocking effects 
by a cold soot layer adjacent to the wall. 
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Task 4. Measurement of Charring Fuel Smoke Points (J. de Ris, 
Xiaofang Cheng) 


The objective of this task is to measure the smoke points of 
several charring solid materials by moving horizontal samples under 
a controlled CO, laser heat source, and to compare results with 
smoke points obtained by the existing vial technique, and correlate 
both results against measurements of incompleteness of combustion. 

It is now well established that the smoke point of a material 
provides a direct measure of a material’s flame sootiness and it is 
anticipated that the smoke point may also provide various other 
non-thermal damage! characteristics of a material. The smoke point 
also correlates the radiant fraction of the total heat release rate 
and the incompleteness of combustion of turbulent diffusion flames 
typical of fires. Clearly, the smoke point is a very important 
material flammability property which needs to be measured for 
practical fuels. It can readily be measured for gaseous, liquid 
and non-charring solid fuels which can be made to support a con- 
trolled steady laminar “candle-like" diffusion flame. The smoke 
point is defined as the maximum flame height which just does not 
emit smoke from the flame tip. A typical 10 cm high smoke-point 
flame requires a fuel supply rate of only 4x10™ g/s. The small 
amount of material required for such measurements should make the 
smoke-—point measurement particularly attractive to chemical indus- 
try researchers involved in developing new materials. 

The major problem with using the smoke point as a material 
flammability property is the difficulty in measuring it for char- 
ring fuels. The pyrolysis rates of charring fuels are usually 
transient. The char produced during transient pyrolysis may 
chemically react with the pyrolysis gases causing them to change 
chemically and alter their smoke-point value. There is also the 
possibility that the pyrolysis vapors may change chemically when 


coming into contact with a heated surface before entering the 


‘Damage due to CO, soot particles and corrosive chemical compounds. 
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laminar candle-like flame. These issues need to be resolved as 
part of the development of a practical smoke-point test method. 

Previously, as part of this grant, M.A. Delichatsios estab-— 
lished that it is possible to measure the smoke-point of a charring 
fuel. He placed samples of PMMA, wood and neoprene in a small 
quartz vial and exposed the vial to a controlled quartz—halogen 
radiant heat source in our flammability research laboratory. By 
Manually increasing and decreasing the radiant exposure, he ob- 
tained a preliminary smoke-point measurement for all three fuels. 
The PMMA value agreed with previous results obtained by Tewarson. 

We have recently assembled an apparatus for measuring the 
smoke point of a charring material by slowly moving a horizontal 
slab under a focused beam from a CO, laser. By adjusting either 
the horizontal feed rate or CO, laser power we can control the 
height of the steady laminar flame. The flame height is measured 
by video camera while the smoke release is measured by its extinc-— 
tion of an electrically chopped IRLED optical beam passing across 
an exhaust chimney above the flame. A threshold circuit powers a 
light bulb seen by the video camera when the smoke release exceeds 
some specified small value. 

Initial measurements of the smoke point for PMMA and two types 
of particle board show that the smoke-point flame height is inde- 
pendent of the horizontal feed rate. Greater feed rates require 
more laser power. A laser power of less than 70 W should be 
sufficient for all materials of interest. 

The laser technique has the particular advantage of allowing 
the fuel gases to immediately enter the attached laminar flame 
without coming into contact with any neighboring hot surface, such 
as the mouth of a heated vial or other fuel holder. The flame 
remains well anchored to the fuel surface because of the intense 
flame to solid heat transfer characteristic of creeping flames. 

We have now tested a new apparatus that can fully exploit this 
laser technique. The apparatus, shown in Figure 1, consists of a 
Stationary glass-—enclosure with perforated metal sheets on two 


sides and the top for air access. Experiments with fuel-gas 
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instead of a solid sample have proven the feasibility of stabiliz-— 
ing the laminar flame with parallel, perforated horizontal plates 
near the flame base, as suggested by a new theory for pulsating 
pool fires. Soot production when the flame height is above the 
smoke point is detected by a near-infrared LED beam while flame 
height is measured with a video camera. 

Some results for PMMA and particleboard are shown in Figures 


2, 3 and 4. A report is in preparation. 
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